A D-AU*  5%5  AIR  FORCE  ROCKET  PROPULSION  LAB  EDWARDS  AFB  CA  F/0  21/9.2 

AN  ANALYSIS  OF  SOLID  PROPELLANT  TRANSIENT  VISCOELASTIC  RESPONSE— ETC <U> 
JAN  82  D  I  THRASHER,  E  CORBETT 

UNCLASSIFIED  AFRPL-TR-81-55  NL 


FILE  COPY 


AFRPL-TR-81-55 


% 


AN  ANALYSIS 
VISCOELASTIC 


OF  A  SOLID  PROPELLANT 
RESPONSE  UNDER  MOTOR 


TRANSIENT 

IGNITION  CONDITIONS 


fi 

Authors:  Durwood  I.  Thrasher 

^-4  Eileen  Corbett 


January  1982 


Final 


report 


for  the  period  9  July 


1980  to  15  April  1981 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED 

The  AFRPL  Technical  Services  Office  has  reviewed  this  report, 
and  it  is  releasable  to  the  National  Technical  Information 
Service,  where  it  will  be  available  to  the  general  public, 
including  foreign  nationals. 


Prepared  for 

AIR  FORCE  ROCKET  PROPULSION  LABORATORY 
DIRECTOR  OF  SCIENCE  AND  TECHNOLOGY 
AIR  FORCE  SYSTEMS  COMMAND 
EDWARDS  AFB  ,  CALIFORNIA  93523 


DT1C 


A 


/ 


l 


NOTICES 

When  U.S.  Government  drawings,  specifications,  or  other  data  are  used  for  any  purpose 
other  than  a  definitely  related  Government  procurement  operation,  the  fact  that  the 
Government  may  have  formulated,  furnished,  or  in  any  way  supplied  the  said  drawings, 
specifications,  or  other  data,  is  not  to  be  regarded  by  implication  or  otherwise,  or  in  any 
manner  licensing  the  holder  or  any  other  person  or  corporation,  or  conveying  any  rights  or 
permission  to  manufacture,  use,  or  sell  any  patented  invention  that  may  be  related  thereto. 

FOREWORD 


This  report  describes  the  work  conducted  by  the  Mechanical  Behavior  and  Aging  Section 
(MKPB)  of  the  Air  Force  Rocket  Propulsion  Laboratory,  Edwards  AFB,  California  93523 
under  Job  Order  No.  2307 M I MD. 

This  Technical  Report  is  approved  for  release  and  distribution  in  accordance  with  the 
distribution  statement  on  the  cover  and  on  the  DD  Form  1473. 

R.  dpHN  MOSS,  Capt,  USAF 
Chief,  Mechanical  Behavior  and 
Aging  Section 


DURWOOD  I.  THRASHER 
Project  Manager 


FOR  THE  DIRECTOR 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (TWl«n  Omit  Bntarad) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS  j 

^BEFORE  COMPLETING  FORM  1 

L  REPORT  NUMBER  2.  GOVT  ACCESSION  NO. 

AFRPL-TR-81-55  ,4f)'/4//  V  ?( 

jfy  |  anm  |  Mirrm— 

«.  TITLE  (and  Submit) 

An  Analysis  of  Solid  Propellant  Transient  Visco¬ 
elastic  Response  Under  Motor  Ignition  Conditions 

5.  TYPE  OF  REPORT  t  PERIOD  COVERED 

Final  Report 

9  Jul  80  to  15  Apr  81 

6.  PERFORMING  ORG.  REPORT  NUMBER 

7.  AUTHORfA) 

Durwood  I.  Thrasher 

Eileen  Corbett 

6.  CONTRACT  OR  GRANT  NUMBERf*> 

9  PERFORMING  ORGANIZATION  NAME  AND  AOORe'sS 

Air  Force  Rocket  Propulsion  Laboratory  (MKPB) 
Edwards  AFB  CA  93523 

10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  &  WORK  UNIT  NUMBERS 

Project  573013  MD 

II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

12.  REPORT  DATE 

January  1982 

13  NUMBER  OF  PAGES 

33 

!«•  MONITORING  AGENCY  NAME  &  AOORESSf/f  dlt/erent  from  Controlling  Olltca) 

IS.  SECURITY  CLASS,  (ol  thla  r abort) 

UNCLASSIFIED 

15a.  DECL  ASSI  FI  CAT!  OnVdOWN  GRADING 

SCHEDULE 

16  Distribution  statement  (oi  thla  Report) 


Approved  for  Public  Release;  Distribution  Unlimited 


\  *7  DISTRIBUTION  STATEMENT  ( of  the  abstract  entered  In  Block  20,  II  different  from  Report) 


i ®  Supplementary  notes 


19  KEY  WORDS  (Continue  on  reverse  side  If  necessary  and  Identify  by  block  number) 

Solid  Rocket  Propellants  Motor  Ignition  Pressurization 

Solid  Propellant  Mechanical  Behavior  Motor  Structural  Analysis 

Sinear  Viscoelasticity 
Transient  Loading 

20.  ABSTRACT  (Continue  on  reverse  elde  If  necessary  and  Identify  by  block  number) 

Programmable  calculator  programs  were  written  to  calculate  stress  responses  for 
stepwise  and  rampwise  approximations  to  a  given  strain  history.  These  pro¬ 
grams  were  tailored  to  calculate  the  error  relative  to  the  exact-history  result 
and  also  to  investigate  different  techniques  of  applying  stepwise  approxima¬ 
tions  and  different  techniques  for  selecting  stress  output  times  within  each 
step  or  ramp.  The  optimum  rampwise  approximation  produced  one-tenth  the  rms 
error  in  stress  produced  by  the  optimum  stepwise  approximation  for  the  same 


DO  ,  ;2nM71  1473  EDITION  OF  I  NOV  66  IS  OBSOLETE  UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  fWlen  Data  Bntarad) 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAQEfWTi^i  DM,  gnWftf;  _ 

number  of  time  Intervals.  Using  only  10  time  Interval's  to  approximate  a  fourth 
order  polynomial  strain  history,  the  rms  error  In  stress  with  the  rampwlse 
approximation  was  only  0.5%.  The  ramowlse  approximation  program  was  then 
rewritten  for  analysis  of  an  arbitrary  Input  strain  history  (up  to  35  ramps). 

A  dynamic  pressurization-strain  test  reported  in  AFRPL-TR-78-68  was  analyzed 
using  the  rampwlse  strain  history  approximation  as  well  as  a  polynomial 
approximation.  For  both  strain  history  approximations,  the  stress  history 
calculated  using  relaxation  data  agreed  well  with  the  measured  stress  history 
up  to  the  peak  stress  point.  The  error  at  peak  stress  for  both  approximations 
was  less  than  20  psl,  or  roughly  2.5%  of  the  measured  value  of  810  psl.  This 
Is  considered  excellent  agreement  with  the  experimental  data. 


UNCLASSIFIED 

security  classification  of  this  faoe dm* 


TABLE  OF  CONTENTS 


Section 

1.  INTRODUCTION 

2.  SUMMARY 

3.  ANALYSIS  DETAILS 

3.1  Overall  Approach 

3.1.1  Power  Law  Modulus  Program 

3.1.2  Polynomial  Program 

3.1.3  Stepwise  Approximation  Program 

3.1.4  Piecewise-Linear  Approximation  Program 

3.1.5  Application  of  Programs 

3.2  Calculations  for  Ramp  Strain  History 

3.3  Calculations  for  Polynomial  Strain  History 

3.4  Application  to  Actual  Data 

4.  CONCLUSIONS  AND  RECOMMENDATIONS 


1 


LIST  OF  FIGURES 


Figure 

Page 

1. 

Step  Approximation  to  Ramp  Strain  History 

11 

2. 

Ramp  Approximation  to  Polynomial  Strain  History 

13 

3. 

Ramp  Strain  History  and  "Exact"  Stress  Response 

15 

4. 

Variation  of  Stress  Error  with  Ko  for  Ke  =  0 
(Step  Approximation  to  Ramp,  30  Time  Intervals) 

17 

5. 

Variation  of  Stress  Error  with  Ka  for  Ke  =  0.5 
(Step  Approximation  to  Ramp,  30  Time  Intervals) 

17 

6. 

Variation  of  Stress  Error  with  Ko  for  Ke  =  1.0 
(Step  Approximation  to  Ramp,  30  Time  Intervals) 

17 

7. 

Variation  of  Stress  Error  with  Number  of  Time  Intervals 
for  Ko  =  0.5  (Step  Approximation  to  Ramp) 

17 

8. 

Variation  of  Stress  Error  with  Number  of  Time  Intervals 
for  Ko  =  0.35  (Step  Approximation  to  Linear  Ramp) 

18 

9. 

Variation  of  Stress  Error  with  Strain  History  Error 
for  Ko  =  0.5  (Step  Approximation  to  Ramp) 

19 

10. 

Variation  of  Stress  Error  with  Strain  History  Error 
for  Optimum  Ko  (Step  Approximation  to  Ramp) 

19 

11. 

Polynomial  Strain  History  Used  in  Study 

20 

12. 

Variation  of  Stress  Error  with  Ko  (Ramp  Approximation 
to  Polynomial,  30  Time  Intervals) 

21 

13. 

Variation  of  Stress  Error  with  Number  of  Time  Intervals 
(Ramp  and  Step  Approximation  to  Polynomial) 

23 

14. 

Strip  Chart  Data  for  Run  #1 

25 

15. 

Pressure-Time  History  Scaled  from  Strip  Chart 

26 

16. 

Strain  and  Stress  Data  Scaled  from  Strip  Chart 

27 

17. 

Relaxation  Modulus  from  ASALM  Uniaxial  Tests  at  2.5% 
Strain  (Modified  Power  Law  and  Log-Quadratic  Approxi¬ 
mations) 

28 

2 


LIST  OF  FIGURES  (concluded) 


Figure  Page 

18.  Relaxation  Modulus  from  ASALM  Uniaxial  Tests  at  2.5/5 

Strain  (Modified  Power  Law  Number  1)  28 

19.  Relaxation  Modulus  from  ASALM  Uniaxial  Tests  at  2.5% 

Strain  (Modified  Power  Law  Number  2)  30 

20.  Polynomial  Strain  History  Results  31 

21.  Piecewise-Linear  Strain  History  Results  31 


3 


LIST  OF  TABLES 

Table  Page 

1  Matrix  of  Parameters  for  Initial  "Ramp"  Calculations  16 

2  Matrix  of  Initial  Ko  Values  Used  In  Piecewise-Linear 

Analysis  22 

3  Matrix  of  Additional  Analysis  Parameters  22 


- 


INTRODUCTION 


The  objective  of  this  project  was  to  evaluate  two  alternative  approaches 
(step  approximation  of  strain  history  and  ramp  approximation  of  strain  history)  to 
calculating  the  transient  stress  response  of  a  linear  viscoelastic  material  to  an 
arbitrary  strain-time  input,  and  to  determine  the  better  approach  for  calculating 
effective  modulus  for  motor  ignition  pressurization  from  stress  relaxation  data. 

The  project  began  with  the  plan  of  expanding  the  capability  of  an  existing 
programmable  calculator  program  which  calculates  stress  response  to  piecewise- 
I inear  strain  histories  so  that  it  could  handle  at  least  20  individual  constont-strain- 
rate  segments.  This  proaram  would  then  be  used  to  calculate  stress  response  using 
both  a  series  of  step  functions  and  a  series  of  ramp  functions  to  approximate  a 
specific  "exact"  continuous-slope  transient  strain  history. 

It  was  found  that  writing  new  calculator  programs  was  more  effective  than 
expanding  the  old  program.  One  of  these  new  programs  can  handle  a  piecewise- 
linear  strain  history  consisting  of  up  to  50  constant-strain-rate  segments  (however, 
only  35  ramps  can  be  loaded  from  the  tape  data  file  in  the  present  data  format). 
Other  calculator  programs  (discussed  later  in  this  report)  were  written  to  perform 
specific  calculations. 

A  study  was  marie  to  refine  the  methods  of  approximating  the  strain  history 
in  each  approach.  To  determine  the  accuracy  of  the  stress  response,  the  stress 
responses  were  compared  with  the  exact  solutior.  for  a  specific  transient  strain 
history. 

2.  SUMMARY 

The  modified  power  law  representation  of  relaxation  modulus  is  amenable  to 
straightforward  calculation  of  linear  viscoelastic  response  to  a  constant  rate 
framp)  strain  input.  A  piecewise-linear  function  can  approximate  a  complex  strain 
history  with  reasonable  accuracy.  This  piecewise-linear  function  is  merely  the  sum 
of  a  series  of  ramp  functions:  the  response  of  the  material  to  the 


piecewise- 1  inear  strain  history  is,  in  turn,  sirrply  the  summation  of  the  responses  to 
the  individual  ramps. 

This  concept  was  applied  successfully  in  the  analyses  described  in  this  report. 
The  rms  error  in  calculated  stress  using  as  few  as  10  ramps  to  approximate  a 
fourth-order-polynomial  strain  history  (quite  similar  to  the  strain  history  typically 
produced  by  motor  iqnition  pressurization)  is  shown  to  be  less  than  0.5%.  This 
accuracy  is  superior  to  that  ohtained  with  100  steps  in  the  usual  stepwise 
approximation. 

The  technique  was  also  applied  to  actual  propellant  data  from  low 
temperature,  hiqh  rate  transient  tests  in  which  the  propellant  strain  was  driven  by 
a  pressure  dynamically  applied  to  a  fluid  surrounding  the  test  specimens.  Using 
stress  relaxation  data  the  calculated  stress  agrees  with  the  measured  stress  within 
?.5%  at  peak  stress.  This  excellent  aqreement  provides  evidence  that  linear 
viscoelasticity  may  be  successfully  applied  to  the  iqnition  pressurization  problem, 
although  the  problem  of  prestrain  effects  was  not  addressed.  Also,  the  propellant 
behavior  departs  dramatically  from  linear  viscoelasticity  once  the  peak  stress  is 
passed;  the  linear  prediction  shows  a  much  hiqher  stress  than  is  actually  observed. 
This  departure  from  the  linear  prediction  may  not  be  important  in  predicting 
failure,  however. 

3.  ANALYSIS  OFTAILS 

The  analysis  was  done  using  a  Hewlett-Packard  9RI5A  programmable 
calculator.  During  this  study  several  programs,  discussed  later,  were  devised  to 
perform  different  calculations  during  the  study.  Two  methods  of  approximating  a 
smooth  strain  time  history  were  compared:  approximation  by  a  series  of  steps  and 
approximation  by  a  series  of  ramps.  The  smooth  curve  used  was  a  fourth-order 
polynomial. 
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3.1  Overall  Approach 

The  calculator  programs  developed  and  their  use  ore  described  below. 

3.1.1  Power  Law  Modulus  Program  -  The  Power  Low  Modulus  Program  was 
written  to  calculate  constants  for  a  modified  power  law  modulus  equation 
applicable  to  an  isothermal  loading  situation,  using  given  relaxation  modulus  data 

as  a  function  of  fime.  The  program  forces  the  power  law  to  fit  three  points  [£.,  f  ({, .)] 
where  5  is  the  "reduced  time"  t/ay,  and  Er  is  relaxation  modulus. 

The  power  law  equation  is 
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Equations  3,  4  and  5  are  solved  in  the  calculator  program  by  iterating  on  an 


assumed  value  of  E^.  The  value  of  E„  converges  quickly  to  within  I  x  l(H>  psi. 

3.1.2  "Polynomial"  Program  -  In  this  program  the  "hereditary  integral"  is 
evaluated  for  a  polynomial  strain  history  given  a  modified  power  law 
representation  of  the  linear  viscoelastic  relaxation  modulus.  The  program  then 
calculates  the  strain  history  and  stress  response. 

The  program  first  determines  the  constants  in  the  polynomial 
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where  T  is  the  time  variable  describing  the  strain  history.  The  constants  are 
determined  from  user  input  values  of  (t,e)  at  the  point  of  maximum  strain  rate  and 
(t,  e  )  at  the  point  of  maximum  strain.  The  stress  at  time  t  is  given  by  the 
hereditary  integral: 
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Evaluation  of  this  integral  using  Equation  6  for  the  strain  history  results  in 


(8) 
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Equations  6,  8  and  9  are  implemented  in  subroutines  which  are  also  used  in  the 
programs  discussed  iri  Sections  3.1.3  and  3.1.4. 

The  calculator  program  pi  ints  out  time,  strain,  and  stress  at  time  defined  by 
a  user-input  maximum  time  and  number  of  increments. 

3.1.3  Stepwise  Approximation  Program.  -  This  program  calculates  a  stepwise 
approximation  to  the  polynomial  strain  history,  evaluates  the  hereditary  integral  to 
obtain  the  stress  response,  and  calculates  root-mean-square  (rms)  errors  for  the 
strain  history  and  stress  response  as  compared  to  those  (the  "exact"  values)  for  the 
polynomial.  Subroutines  from  the  "Polynomial"  program  are  used  to  calculate  the 
"exact"  stress  and  strain  values. 

The  power  law  used  is  Equation  2.  The  hereditary  integral,  equation  7,  when 
evaluated  at  time  tj  during  the  step,  becomes 
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and  where  (ep)f  refers  to  the  strain  given  by  the  polynomial  (Eq  6)  at  time  t;  *1  is 
the  time  at  which  the  jth  step  occurs;  and  the  parameters  K  and  K0  are  defined 
pictorially  in  Figure  I.  (While  Figure  I  shows  the  actual  strain  history  as  a  ramp 
function,  the  definition  of  and  Kq  are  independent  of  the  actual  strain 

history.)  Variations  in  and  Kq  allow  different  schemes  for  defining  the 
stepwise  strain  history  (given  a  "real"  continuous  history)  and  computing  a 
corresponding  stress  history. 

The  program  implements  Equations  10,  II,  and  12,  calculates  the  rms  relative 
error  in  strain  at  m  evenly  spaced  points  in  the  strain  history  (relative  to  the 
polynomial  history)  and  calculates  the  rms  relative  error  in  the  step-response 
stresses  relative  to  the  "exact"  polynomial-history  response  (only  the  stress  values 
at  times  tj  are  considered).  The  user  controls  the  number  of  steps  (NSTEPS)  and 
the  parameter  m  as  well  as  K£  and  K^. 

3.1.4  Piecewise-linear  Approximation  Program  -  This  program  calculates  a 
piecewise-linear  approximation  to  the  polynomial  strain  history,  evaluates  the 
hereditary  integral,  and  calculates  rms  relative  errors  for  the  strain  history  and 
stress  response.  Subroutines  from  the  "Polynomial"  program  are  implemented  to 
calculated  exact  stress  and  strain  values. 

The  power  law  used  is  Equation  2.  The  hereditary  integral,  Equation  7,  when 

evaluated  at  time  tj  during  the  jth  step,  becomes: 
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\ Actual  strain  history 
( Ramp :  i  =  Rt ) 


Stupwist  approximation  to  strain  history 
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where  («p)t  refers  to  the  strain  given  by  the  polynomial  (Equation  6)  at  time  t,  j  is 
the  time  at  which  the  jth  ramp  begins,  and  the  parameter  is  defined  in  Figure  2. 


The  program  implements  Equations  13,  14  and  15.  The  strain  history  and 
stress  response  errors  are  calculated  the  same  way  as  in  the  stepwise 
approximation  program  (Section  3.1.3).  The  user  inputs  the  number  of  steps 
(NSTEPS),  the  parameter  m,  and 

3.1.5  Application  of  Programs  -  The  programs  described  in  Sections  3.1.3  and 
3.1.4  were  used  first  to  find  the  optimum  values  of  Kt  and  Ka  for  approximating  a 
ramp  (taken  as  apolynomial  with  zero  constants  except  for  the  linear  term)  with 
steps,  and  then  to  find  the  optimum  value  of  K for  approximating  the  polynomial 
strain  history  with  ramps.  Finally,  the  performance  of  stepwise  approximations  to 
the  polynomial  and  piecewise- 1  inear  approximations  to  the  polynomial  were 
compared  to  select  the  best  way  of  approximating  an  arbitrary  history.  The 
piecewise-l inear  approximation  was  found  to  be  best.  These  calculations  are 
discussed  in  detail  in  Sections  3.2  and  3.3. 

Finally,  a  new  progrcm,  "Ramp-Series  Stress  Response  for  Arbitrary  Load 
History,"  was  written  to  calculate  a  piecewise-linear  approximation  to  an  arbitrary 
strain  history  and  to  the  corresponding  stress  response.  The  new  program  is  very 
similar  to  the  "Piecewise-linear  Approximation  Program."  The  primary  differences 
are  that  the  new  program  uses  a  discrete  table  of  time  and  strain  values  instead  of 
the  polynomial  history,  and  that  it  does  not  contain  the  error  calculations.  This 
new  program  was  applied  to  an  actual  strain  history  from  data  reported  in  AFRPL- 
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TR-78-68,  "Improved  Solid  Propellant  Mechanical  Properties  Measurement  for 
Structural  Analysis  Input."  The  resulting  stress  history,  as  well  as  the  stress 
history  resulting  from  a  polynomial  strain  history  which  closely  matches  the  actual 
strain  history,  were  compared  with  the  actual  measured  stress  history.  These 
results  are  discussed  in  Section  3.4. 

3.2  Calculations  for  Ramp  Strain  History 

An  arbitrary  modulus  law  (Equation  2)  was  used  with  Ec  =  2.3859  x  l(P 
psi  (16.450  MPa),  E  =  55.408  psi  (0.38202  MPa),  and  n  =  -0.32790,  with  t  in 

OO 

minutes  (for  t  in  seconds,  the  value  of  E0  would  be  623.16  psi  (4.297  MPa)).  The 
Stepwise  Approximation  Program  was  used  with  a  ramp  substituted  for  the 
polynomial  (C  =  21.923,  C2  =  C3  =  C4  =  0,  for  t  in  minutes;  C|  =  0.36538,  C2  =  C3  = 
C4  =  0  for  t  in  seconds). 

Figure  3  shows  the  ramp  strain  history  and  the  resulting  "exact"  stress 
response  used  in  this  part  of  the  study. 

Initial  runs  with  the  Stepwise  Approximation  Program  were  made  using 

selected  values  of  the  parameters  K  e,  «a,  and  NSTEPS  as  shown  in  Table  I.  The 

results  (shown  in  Figures  4-6)  show  the  rms  stress  error  as  a  function  of  K  for  K  = 

0  e 

0,  0.5,  and  1.0  (these  Ke  values  were  chosen  because  they  are  simple  to  use  for 

evenly  spaced  time  intervals  when  the  data  is  available  only  in  discrete  form). 

Additional  runs  were  made  to  search  out  the  optimum  (i.e.,  least-error)  value  of  K 

0 

for  each  K  .  The  best  results  were  obtained  with  K  =  0.5  and  K  =  0.35. 
t  e  e 

Interestingly,  it  can  be  shown  that,  for  a  single  step  approximating  a  ramp  over  the 
time  interval  Tm,  the  value  of  which  yields  identical  stresses  for  the  ramp  and 
the  step  (for  Ke  0)  is 


K  =  K 

a  t 


r*  *  n 
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Figure  3. 
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Ramp  Strain  History  and  "Exact"  Stress  Response. 
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TABLE  1.  MATRIX  OF  PARAMETERS  FOR  INITIAL  "RAMPS"  CALCULATIONS 


](> 


Ka  (dimensionless) 


K,  (dimensionless) 


Figure  4.  Variation  of  Stress  Error 
with  K0  for  K£  =  0  (Step  Approximation 
to  Ramp,  30  Time  Intervals). 


K  *  (dimensionless) 

Figure  6.  Variation  of  Stress  Error 
with  Kq  for  Kj.  =  1.0  (Step  Approxima¬ 
tion  to  Ramp,  30  Time  Intervals) . 


Figure  5.  Variation  of  Stress  Error 
with  K0  for  K?  =  0.5  (Step  Approximation 
to  Ramp,  30  Time  Intervals). 


Figure  7.  Variation  of  Stress  Error 
with  Number  of  Time  Intervals  for 
Ka  =  0.5  (Step  Approximation  to  Ramp) 


For  t  m  -  1.3  x  I  O'3  and  the  'modulus  law  constants  used  in  the  study,  this  equation 
yields  K*  =  0.674  K(  ,  which  agrees  fairly  closely  with  the  results  in  Figures  5  and 
6.  The  optimum  value  of  K0  may  therefore  be  a  function  of  the  modulus  constants. 
This  factor  was  not  explored  further  in  the  study. 


Figure  8.  Variation  of  Stress  Error  with  Nunoer  of  Time  Intervals 
for  Kg-  =  .35  (Step  Approximation  to  Linear  Ramp). 

Figures  7  and  8  show  the  variation  of  stress  error  with  NSTEPS.  As  can  be 
verified  by  a  log-log  plot  of  the  variables,  the  stress  error  varies  as  l/_'v  NSTEPS,  so 
that  the  user  pays  a  heavy  price  in  computation  time  to  improve  accuracy  by 
increasing  the  number  of  steps.  It  is  obvious  from  these  figures  that  using  either  Ka 
=  0.5  or  K(  -  1.0  is  vastly  preferable  to  using  K(  =  0,  and  that  using  the  optimum 
value  of  K  for  each  K(  provides  a  substantial  improvement  in  the  stress  error  in 
comparison  to  an  arbitrary  value  of  Ka  =  0.5. 
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Figure  9.  Variation  of  Stress  Error  with  Strain  History  Error  for  K  =0.5 
(Step  Approximation  to  Ramp) . 


Figure  10.  Variation  of  Stress  Error  with  Strain  History  Error  for  Optimum 
(Step  Approximation  to  Ramp) . 

Figures  9  and  10  show  the  variation  of  stress  error  with  strain  history  error; 
the  improvement  provided  by  the  optimum  values  is  again  clearly  evident. 

The  lowest  error  is  obtained  using  K(  =  0.5  and  the  corresponding  value  of  Kff 
(K^  -  0.35).  These  values  were  selected  for  use  in  the  continued  study  described  in 
the  following  section. 
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Figure  11.  Polynomial  Strain  History  Used  in  Study. 


3.3  Calculations  For  Polynomial  Strain  History 

The  arbitrary  modulus  law  (Eq  2)  used  in  the  calculations  for  a  ramp  strain 

A 

history  and  the  constants  E0,  E  .  and  n  are  given  in  Section  3.2  were  also  used  in 
this  part  of  the  study. 

The  Piecewise-linear  Approximation  Program  was  used  with  a  polynomial 

strain  history.  The  constants  were  determined,  using  the  control  points  shown  in 

Figure  I  la,  with  the  "Polynomial"  program.  The  resulting  constants  are: 

(Time  in  Minutes)  (Time  in  Seconds) 

C  |  —  2.7549 

C2  =  6.3493  x  10** 

C3  -  -5.0683  x  I07 

Cz,  ::  1.0142  X  I010 

The  resulting  "exact"  stress  response  is  shown  in  Figure  I  lb. 


C|  =  1.6529  x  10' 
C2  =  2.2857  x  108 
C3  -  -1.0948  x  I013 
-  1.3144  x  10)7 


r\ 


Stress  error  (dimensionless) 


Figure 


K  (j- 


12.  Variation  of  Stress  Error  with  K 
to  Polynomial,  30  Time  Intervals 
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Initial  runs  with  the  Piecewise-linear  Program  were  made  using  selected 
values  of  the  parameters  Kg  and  NSTEPS  as  shown  in  Table  2.  The  results  are 
shown  in  Figure  12. 

TABLE  2.  MATRIX  OF  INITIAL  Ka  VALUES  USED  IN  PIECEWISE-LINEAR  ANALYSES 


.001 

.25 

.5 

.75 

1.0 

NSteps  = 

X 

X 

X 

X 

X 

30 

Based  on  the  discussion  in  Section  3.2,  the  optimum  value  of  Ka  is  probably 
dependent  on  the  modulus  constants.  This  study  did  not  explore  this  factor 
further. 


Figure  13  (based  on  the  results  of  the  additional  runs  indicated  in  Table  3) 
shows  the  piecewise-linear  approximation  is  superior  to  the  stepwise 
approximation.  The  piecewise-linear  approximation  has  a  much  lower  stress  error 


for  any  combination  of  NSTEPS  and  the  controlling  parameter  K^which  forces  the 
stepwise  program's  error  below  5%.  For  example,  using  optimum  values  of  Ka,  the 
piecewise-linear  approximation  ,  ith  NSTEPS  -  10  yields  a  lower  error  than 


Figure  13. 


Variation  of  Stress  Error  with  Number  of  Time 
Intervals  (Ramp  and  Step  Approximations  to  Polynomial) 


usinq  the  stepwise  approximation  with  NSTFPS  -  100.  If  an  rms  error  of  0.5%  is 
acceptable  in  the  stress  response,  then  a  value  of  NSTEPS  =  10  should  be  sufficient 
in  a  piecewise-linear  approximation  of  a  general  strain  history. 

3.4  Application  to  Actual  Data 

AFRPL-TR-78-68,  "Improved  Solid  Propellant  Mechanical  Properties 
Measurement  for  Structural  Analysis  Input"  (Thiokol/Huntsville,  September  1978), 
reports  the  results  of  dynamic  tests  of  propellant  samples  under  simulated  motor 
ignition  conditions.  In  these  tests  propellant  samples  were  immersed  in  a 
temperature-conditioned  fluid  and  held  at  a  predetermined  strain  level  to  simulate 
motor  storage  thermal  loads.  The  fluid  was  then  pressurized  with  gas,  providing  a 
transient  pressure  history  similar  to  a  motor  ignition  pressurization.  The  test 
equipment  was  arranged  so  that  the  pressurized  fluid  drove  the  test  machine  ram 
downward  against  the  resistance  of  steel  springs  (as  well  as  the  propellant  sample's 
resistance  and  the  inertia  of  moving  parts).  The  result  was  a  transient  strain 
history  that  simulated  the  strain  at  the  inner  bore  of  a  rocket  motor  undergoing 
ignition  pressurization. 

Since  the  stress-strain  data  for  some  of  these  tests  was  reported,  along  with 
conventional  relaxation  test  data,  an  ideal  test  case  was  provided  for  the 
techniques  studied  in  the  present  effort.  Fiqure  14  is  a  photograph  of  the 
oscillograph  ("strip  chart")  data  (Figure  9  in  AFRPL-TR-78-68).  This  data 
corresponds  to  Run  No  I  in  Table  F-5  of  AFRPL-TR-78-68.  The  prestrain  level  for 
the  tests  was  zero.  Other  test  parameters  and  results  given  in  Table  F-5  of 
AFRPL-TR-78-68  were: 

Temperature:  -65°F  (-53.9°C) 

Maximum  Pressure:  1310  psi  (9.03  MPa) 

Time  to  Maximum  Pressure;  67.9  ms 

Strain  at  Maximum  Pressui •  ■  0.0317  in/in  (.0317  mm/mm) 

732  psi  (5.05  MPa) 


Stress  at  Maximum  Pressure: 


These-  values  were  used  to  scale  the  data  shown  in  F  igure  14.  In  scaling  the 


displacement  traces  the  amplitude  of  both  !)|  and  F>2  were  taken  as  0.0317  in/in 
(0.0317  rnm/imn);  the  resulting  scaled  strains  were  then  averaged.  The  scaled  data 
is  given  in  figures  15  and  16.  (A  careful  examination  of  F  igure  14  will  reveal  that 

the  image  is  distorted,  presumably  by  perspective  effects  in  the  photography.  An 
effort  was  made  to  accurately  account  for  the  distortion  in  scaling  the  data.) 

Comparison  of  Figure  15  with  Figure  16(a)  shows  that  the  strain  history  lags 
the  pressure  history  somewhat  but  otherwise  corresponds  fairly  closely.  The  strain 
history  has  the  same  general  appearance  as  the  polynomial  history  shown  in  Figure 
I  I  (up  to  the  peak  value),  except  for  an  added  oscillation  or  "bumpiness",  which 
appears  even  more  strongly  in  the  stress  history  (F  igure  16(b)). 
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Figure  16.  Strain  and  Stress  Data  Scaled  from  Strip  Chart. 
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The  relaxation  modulus  for  this  propellant  was  determined  at  several 
different  strain  levels.  The  data  reported  for  2.5%  strain  as  well  as  two  different 
analytical  representations  is  shown  in  Figure  17.  The  power  law  and  the  log- 
quadratic  equations  are  tangent  at  the  log 1 0  (t/aT ^  =  -7  point  (with  t/aT  in 
minutes).  While  the  combination  of  these  two  functions  fits  the  data  well  over  the 
entire  range  of  data,  the  log-quadratic  function  produces  a  hereditary  integral 
which  must  be  numerically  evaluated  for  any  strain  history  other  than  a  step 
function.  In  practice,  the  power  law  is  the  only  part  of  the  combination  involved  in 
evaulating  propellant  response  under  the  loading  conditions  of  interest  in  this 
report  because  only  reduced  times  less  than  I0"7  minutes  are  involved.  It  was  felt 
that  the  Thiokol  power  law  had  a  higher  slope  than  was  supported  by  the  data,  so 
two  alternative  modified  power  laws  were  determined  using  the  Power  Law 
Modulus  Program  discussed  earlier.  These  alternative  representations  are 
"Modified  Power  Law  Number  I",  shown  in  Figure  18,  and  "Modified  Power  Law 
Number  2",  shown  in  Figure  19.  The  difference  between  these  two  representations 
is  in  the  center  control  point  (the  modulus  equation  was  forced  through  the  three 
points  shown  in  each  figure).  The  modulus  constants  (for  a  temperature  of  -65°F 
(-53.9°C))determined  for  the  modified  power  law  equation  (Eq.  2)  for  the  English 
system  of  units  (modulus  in  psi,  t/aj  in  minutes)  are: 


o,  psi 

psi 

n  (dim) 

Modified  Power  Law  No  1 

2497 

38.8 

-0.28252 

Modified  Power  Law  No  2 

2426 

44.9 

-0.30382 

For  the  SI  system  of  units  (modulus  in  MPa,  t/a j  in  seconds),  the  constants  are: 


rr 

o 

TD 

c* 

E=  ,  MPa 

n  (dim) 

Modified  Power  Law  No  1 

54.74 

0.268 

-0.28252 

Modified  Power  Law  No  2 

58.03 

0.310 

-0.30382 

Thp  two  modified  power  In-/'  representations  /ire  rotated  to  the  relaxation 


dntn  hy  the  WLf  shift  factor  (a  )  equation. 
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which  was  determined  hv  Ihiokol  specifically  for  the  relaxation  data  in  Figures  17, 
18,  and  lQ.  In  this  equation.  !  .  in  K.  The  rale  Motions  in  this  report  use  a  value 
of  298.15  where  the  value  2°8  appear  in  F  gentian  16:  this  corresponds  to  a 
reference  temperature  (for  which  Log  af  "  (T)  of  77°l  or  25°C . 

Examination  of  I  iqures  18  and  1°  shows  that  "Modified  Power  Law  I  lumber  2" 
nqrees  better  with  the  measured  data  at  very  short  times  (t/aj  <  10“^)  min),  while 
Modified  Power  Law  Number  I  agrees  bettor  with  the  data  for  lonqer  reduced 
times.  We  are  concerned  with  reduced  times  shorter  than  10"^  minutes,  so 
"Modified  Pov/er  Law  Number  2"  appears  preferable  to  "Modified  Power  Law 
Number  I."  "Modified  Power  l.aw  Number  2"  also  appears  to  be  inore  consistent 
with  the  relaxation  dal  than  does  the  Thiokol  power  law. 
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Figure  20.  Polynomial  Strain  History  Results. 


0  0.02  0.04  0.06  0.08  0.10 

Time,  sec 


(b)  Stress  Response 

Figure  21.  Plecewlse-L Inear  Strain  History  Results. 


The  strain  histories  considered  are  shown  in  Figures  20(a)  and  21(a).  The  data 
points  plotted  as  circles  were  scaled  from  Figure  14  using  strain  at  maximum 
pressure  to  determine  scale  factors. 

The  polynomial  strain  history  was  forced  through  a  point  of  inflection  and  a 
maximum  at  the  control  points  shown  in  Figure  20(a).  The  resulting  constants  in 
Equation  6  were: 

Time  in  Minutes  Time  in  Seconds 

C  |  =  5.5035  x  1 02  (min)- 1  C  |  =  9. 1 725  x  1 0'**  (sec)- 1 

C2  =  8.8909  x  10**  (min)'2  C2  =  24.697  (sec)'2 

C3  =  -7.7605  x  I07  (min)-3  C3  =  -3.5928  x  I02  (sec)'3 

C4  =  1.8332  x  I010  C4  =  1.4145  x  I03  (sec)-** 
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Pie  rosul tine]  stress  ifsponses  tire  shown  in  I  iqure  />()(l>).  along  with  the  stress 
dum  (i  tr<  les)  scaled  from  I  iqure  I  k,  t  <->r  bold  "Modified  Power  low  Number  I"  find 
", Modified  Power  Low  I  lumber  As  expected,  "Modified  Power  Low  Momher  ?" 

o.jrees  Letter  with  the  (Into.  I  lotahly,  the  "bumpiness"  seen  in  the  (Into  is  not 
raptured  in  the  stress  response  (the  polynomial  strain  history  smooths  through  the 
"bumps”  in  r be  strain  tristory).  The  calculated  stress  for  "Modified  Power  Law 
tore, tier  agrees  guite  veil  with  the  measured  stress  through  the  maximum 
measured  tress  point.  However ,  once  the  stress  begins  to  decrease,  the  calculated 
visioelnstii  response  fails  to  drop  oft  ns  quickly  ns  the-  measured  stress.  By  the 
t'  id  peak  pressure  r<‘  !  i  . !  1  x  lit*  •  n  in  I,'."/’  mVc)  the  calculated  stress  is 
11  . i  f  ■n;u: -r  t- ’di  !>'<•■  i  ,t  i iv  r i >. '  ■-  r  r 

e  niecrwise-linr  *rr  s:r  1  1  r •  r  /  wo*  i  >>ri.sfrur  fed  through  the  date;  points 

free  L  ignre  IT,  ns  sf'own  in  I  iqure  Min,).  Ii  >•  value  of  used  (see  I  quatiori  15 
arm  T  iqure  7\  was  b.H)3  (bus  value  v'as  found  to  produce  minimum  root  mean 
sgnere  error  in  the  sin-:/  ]|v  in  Section  >.3.  which  used  a  similar  strain 

history  and  modulus  constants).  I  he  results  were  essenl inllv  identical  to  those  for 
the  polynomial  strain  history,  excop'  ttu;t  !h e  elculoted  stress  rrsponse  foiloweii 
tiif-  "blimps"  seep  in  the  nr  tuul  data. 

P.a-.ed  on  the  results,  a'i-  observe  tt.nl: 

( !  )  "Modified  Power  -  ov  !  lumber  /"  is  a  <|oud  ap|  ro>  imntion  to  the  linear 
visr  oelostit  rnorfulus  at  high  rates  (although  we  bad  In  extrapolate  beyond  the 
range  of  measured  data  fur  the  problem  at  hand). 

(/)  !  inenr  vis'  oelnsticity  (based  on  (tie  limited  bain  available)  appears  to 

ficcuratel  •/  prerlict  stress  response  to  a  transient  strain  history  without  a  prestrain, 
ft! ffroar (ft  stress.-s  ! .evorv!  the  time  of  at  fun!  peak  stress  ore  overpredicted.  (Some 
damage  offer  ts  rnciv  be  the  '  .-use  of  the  more  r  r  if  >i<t  drop-off  in  measured  stress). 

I  h'-  error  at  peak  stress  fur  me  lie  pols'nemiul  and  piecewise  linear  strain  history 
appro  durations  is  less  than  /■)  g-.  1  u(  tl.<  •<  ensured  value  at  MIO  psi,  or  less  than 

t.  . 


Jfi' 


4.  CONCLUSIONS  AND  RECOMMENDATIONS 

This  successful  analysis  effort  fully  met  the  objectives.  Some  innovative 
approaches  to  the  analysis  of  propellant  viscoelastic  response  under  motor  ignition 
pressurization  conditions  were  developed,  and  a  number  of  useful  calculator 

programs  were  created  for  AFRPL  in-house  use. 

We  recommend  that  the  techniques  developed  in  this  program  be  extended  to 
the  more  difficult  problem  of  thermal  transient  loading  with  coupled  strain  and 
temperature  histories.  We  also  recommend  that  the  techniques  be  applied  to 
problems  involving  prestrain  to  determine  whether  prestrain  causes  the  linear 
approach  to  fail. 
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